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Abstract 

In these Lectures we review the Minimal Supersymmetric Standard Model as well 
as some of its extensions that include R-Parity violation. The cases of spontaneous 
breaking of R-Parity as well as that of explicit violation through bilinear terms in 
the superpotential are studied in detail. The signals at LEP and the prospects for 
LHC are discussed. 
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1 The Minimal Supersymmetric Standard Model 



1.1 Introduction and Motivation 

In recent years it has been established with great precision (in some cases better than 
0.1%) that the interactions of the gauge bosons with the fermions are described by the 
Standard Model (SM) ||. However other sectors of the SM have been tested to a much 
lesser degree. In fact only now we are beginning to probe the self-interactions of the 
gauge bosons through their pair production at the Tevatron || and LEP [|j] and the 
Higgs sector, responsible for the symmetry breaking has not yet been tested. 

Despite all its successes, the SM still has many unanswered questions. Among the various 
candidates to Physics Beyond the Standard Model, supersymmetric theories play a special 
role. Although there is not yet direct experimental evidence for supersymmetry (SUSY), 
there are many theoretical arguments indicating that SUSY might be of relevance for 
physics below the 1 TeV scale. 

The most commonly invoked theoretical arguments for SUSY are: 

i. Interrelates matter fields (leptons and quarks) with force fields (gauge and/or Higgs 



ii. As local SUSY implies gravity (supergravity) it could provide a way to unify gravity 
with the other interactions. 

Hi. As SUSY and supergravity have fewer divergences than conventional filed theories, 
the hope is that it could provide a consistent (finite) quantum gravity theory. 

iv. SUSY can help to understand the mass problem, in particular solve the naturalness 
problem ( and in some models even the hierarchy problem) if SUSY particles have 
masses < C(lTeV). 

As it is the last argument that makes SUSY particularly attractive for the experiments 
being done or proposed for the next decade, let us explain the idea in more detail. As 
the SM is not asymptotically free, at some energy scale A, the interactions must become 
strong indicating the existence of new physics. Candidates for this scale are, for instance, 
M x ^ C(10 16 GeV) in GUT's or more fundamentally the Planck scale M P ~ O(10 19 GeV). 
This alone does not indicate that the new physics should be related to SUSY, but the so- 
called mass problem does. The only consistent way to give masses to the gauge bosons and 
fermions is through the Higgs mechanism involving at least one spin zero Higgs boson. 
Although the Higgs boson mass is not fixed by the theory, a value much bigger than 
< H° >~ Gp 1 ^ 2 ~ 250 GeV would imply that the Higgs sector would be strongly coupled 
making it difficult to understand why we are seeing an apparently successful perturbation 
theory at low energies. Now the one loop radiative corrections to the Higgs boson mass 
would give 



bosons) . 




(1) 
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which would be too large if A is identified with Aqut or Api anc k- SUSY cures this problem 
in the following way. If SUSY were exact, radiative corrections to the scalar masses 
squared would be absent because the contribution of fermion loops exactly cancels against 
the boson loops. Therefore if SUSY is broken, as it must, we should have 



6m H = °{^) (2) 

We conclude that SUSY provides a solution for the the naturalness problem if the masses 
of the superpartners are below (9(1 TeV). This is the main reason behind all the phe- 
nomenological interest in SUSY. 

In the following we will give a brief review of the main aspects of the SUSY extension of 
the SM, the so-called Minimal Supersymmetric Standard Model (MSSM). Almost all the 
material is covered in many excellent reviews that exist in the literature ||. 

1.2 SUSY Algebra, Representations and Particle Content 
1.2.1 SUSY Algebra 

The SUSY generators obey the following algebra 

{Qa,Qp} = o 

{Q & ,Qp} = 

{Qa,Qp} = 2(a%P, (3) 

where 

^ = (1,0 ; <? = (l-j) (4) 

and oe,/3,a,$ = 1,2 (Weyl 2-component spinor notation). The commutation relations 
with the generators of the Poincare group are 

[P^Qa] = 

[M^,Q a ] = -i(en a fi Q fi 
From these relations one can easily derive that the two invariants of the Poincare group, 

p2 = PnP a 

(5) 



a 

w 2 = w a w a 

where is the Pauli-Lubanski vector operator 



W, = -^ w ilW CT (6) 
are no longer invariants of the Super Poincare group. In fact 

[Qa,P 2 ] =0 

[Qa,W 2 }^0 (7) 

showing that the irreducible multiplets will have particles of the same mass but different 
spin. 
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1.2.2 Simple Results from the Algebra 

From the supersymmetric algebra one can derive two important results: 

A. Number of Bosons = Number of Fermions 

We have 

Q a \B>=\F> ; (-1) Nf \B >= \B > 
Q a \F>=\B> ; (-1) Nf \F >= -\F > 
where (—1)^ is the fermion number of a given state. Then we obtain 

Q„(-l)"' = -{-l) N 'Q a 
Using this relation we can show that 



(9) 



Tr 



■l) N *{Qa,Q*} 



Tr 
Tr 




-l) N 'Q a Q& + (-l) N *Q & Q 

-Q Q (-lf F Q.+Q a (-lf F Q 



But using Eq. @ we also have 

Tr ((-1)"' {Q ai Q a }) = Tr ((-1)^2^ 
This in turn implies 



(10) 



Tr(— l)jv P = #Bosons — ^Fermions = 



showing that in a given representation the number of degrees of freedom of the 
bosons equals those of the fermions. 

B. (0|#|0) > 

From the algebra we get 

{Qi,Qi} + {Q2,Q2} = 2TrK)P M 

= AH 

Then 



H = \ (QiQi + Q2Q2 + QiQi + Q2Q2 
and 

(0|#|0) = [\\Qx |0) || 2 + \\Qx |0) || 2 + 11$! |0> || 2 + WQi |0> || 2 

> 

showing that the energy of the vacuum state is always positive definite. 



(11) 
(12) 



(13) 



3 



1.2.3 SUSY Representations 

We consider separately the massive and the massless case. 

A. Massive case 

In the rest frame 

{Q«,Qd} =2m <U (14) 

This algebra is similar to the algebra of the spin 1/2 creation and annihilation 
operators. Choose \Q) such that 

Q 1 \Q)=Q 2 \n) = (15) 

Then we have 4 states 

\Q) ; Q x \n) ; Q 2 \Q) ; Q,Q 2 \Q) (16) 

If J3 |Q) = J3 \Q) we show in Table [l] the values of J3 for the 4 states. We notice 



State 


J 3 Eigenvalue 


|n> 




Qiin) 


J3 + I 


g 2 |n) 


J3 - | 


g^in) 


J3 



Table 1: Massive states 

that there two bosons and two fermions and that the states are separated by one 
half unit of spin. 

B. Massless case 

If m = then we can choose P M = (E, 0, 0, E). In this frame 

{Q a ,Q & } = M ad 
where the matrix M takes the form 



M 




AE 



Then 



{Q 2 ,Q 2 }=4E 
all others vanish. We have then just two states 

in) ; Q 2 |n) 

If J 3 |n) = A |n) we have the states shown in Table ||, 



(17) 

(18) 
(19) 
(20) 
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State 


J3 Eigenvalue 


Q 2 |n> 


A 

A-i 



Table 2: Massless states 
1.3 How to Build a SUSY Model 

To construct supersymmetric Lagrangians one normally uses superfield methods (see for 
instance f|). In these lectures we do not have time to go into the details of that construc- 
tion. Therefore we will take a more pragmatic view and give the results in the form of a 
recipe. To simplify matters even further we just consider one gauge group G. Then the 
gauge bosons W° are in the adjoint representation of G and are described by the massless 
gauge supermultiplet 

V a = (A a , Wp (21) 

where A a are the superpartners of the gauge bosons, the so-called gauginos. We also 
consider only one matter chiral superfield 

= ; (i = l,...,N) (22) 

belonging to some N dimensional representation of G. We will give the rules for the 
different parts of the Lagrangian for these superfields. The generalization to the case 
where we have more complicated gauge groups and more matter supermultiplets, like in 
the MSSM, is straightforward. 

1.3.1 Kinetic Terms 

Like in any gauge theory we have 

C Un = -\F; u F a ^ + % - Y^D^ + (D,A) j D^A + i^D.P^ (23) 
where the covariant derivative is 

D, = d, + igW;T a (24) 
In Eq. (|23|) one should note that ip is left handed and that A is a Majorana spinor. 

1.3.2 Self Interactions of the Gauge Multiplet 

For a non Abelian gauge group G we have the usual self-interactions (cubic and quartic) 
of the gauge bosons with themselves. These are well known and we do write them here 
again. But we have a new interaction of the gauge bosons with the gauginos. In two 
component spinor notation it reads || 

Cxxw = Wfabc AVA 6 W; + h.c. (25) 
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where f a b c are the structure constants of the gauge group G and the matrices were 
introduced in Eq. (f|). 

1.3.3 Interactions of the Gauge and Matter Multiplets 

In the usual non Abelian gauge theories we have the interactions of the gauge bosons 
with the fermions and scalars of the theory. In the supersymmetric case we also have 
interactions of the gauginos with the fermions and scalars of the chiral matter multiplet. 
The general form, in two component spinor notation, is 0, 

U w = -gT«w; + iAfipAj) + g 2 (t^)^ W*W» h A* A, 

+zgV2T«(\ a ^A*-\%A J ) (26) 

where the new interactions of the gauginos with the fermions and scalars are given in the 
last term. 



1.3.4 Self Interactions of the Matter Multiplet 



These correspond in non supersymmetric gauge theories both to the Yukawa interactions 
and to the scalar potential. In supersymmetric gauge theories we have less freedom to 
construct these terms. The first step is to construct the superpotential W. This must be 
a gauge invariant polynomial function of the scalar components of the chiral multiplet $j, 
that is Ai. It does not depend on A*. In order to have renormalizable theories, the degree 
of the polynomial must be at most three. This is in contrast with non supersymmetric 
gauge theories where we can construct the scalar potential with a polynomial up to the 
fourth degree. 

Once we have the superpotential W, then the theory is defined and the Yukawa interac- 
tions are 

d 2 W , , ( d 2 W Y -- 



^Yukawa 

and the scalar potential is 
where 



dAtdAj " TJ \dAidAj 



V. 



scalar 

F, 
D a 



1 D a D a + FiF* 



(27) 
(28) 



dW 
dA~ 

g A*T°jAj 



(29) 



We see easily from these equations that, if the polynomial degree of W were higher than 
three, then the scalar potential would be a polynomial of degree higher than four and 
hence non renormalizable. 



6 



1.3.5 Supersymmetry Breaking Lagrangian 

As we have not discovered superpartners of the known particles with the same mass, we 
conclude that SUSY has to be broken. How this done is the least understood sector of 
the theory. In fact, as we shall see, the majority of the unknown parameters come from 
this sector. As we do not want to spoil the good features of SUSY, the form of these 
SUSY breaking terms has to obey some restrictions. It has been shown that the added 
terms can only be mass terms, or have the same form of the superpotential, with arbitrary 
coefficients. These are called soft terms. Therefore, for the model that we are considering, 
the general form would beQ 

C SB = ml Re(A 2 ) + m\ Im{A 2 ) - m 3 (\ a \ a + XT) + m 4 {A 3 + h.c.) (30) 

where A 2 and A 3 are gauge invariant combinations of the scalar fields. From its form, we 
see that it only affects the scalar potential and the masses of the gauginos. The parameters 
rrii have the dimension of a mass and are in general arbitrary. 

1.3.6 R-Parity 

In many models there is a multiplicatively conserved quantum number the called R-parity. 
It is defined as 

R = (_1)2J+3B+L (31) 

With this definition it has the value +1 for the known particles and —1 for their su- 
perpartners. The MSSM it is a model where R-parity is conserved. The conservation 
of R-parity has three important consequences: i) SUSY particles are pair produced, ii) 
SUSY particles decay into SUSY particles and in) The lightest SUSY particle is stable 
(LSP). In Sections 2 and 3 we will discuss models where R-parity is not conserved. 

1.4 The Minimal Supersymmetric Standard Model 
1.4.1 The Gauge Group and Particle Content 

We want to describe the supersymmetric version of the SM. Therefore the gauge group is 
considered to be that of the SM, that is 

G = SU C (3)®SU L (2)®U Y (1) (32) 

We will now describe the minimal particle content. 

• Gauge Fields 

We want to have gauge fields for the gauge group G = SU C (3) (g) SUl{2) ® Uy{l). 
Therefore we will need three vector superfields (or vector supermultiplets) Vi with 
the following components: 

2 We do not consider a term linear in A because we are assuming that $, and hence A, are not gauge 
singlets. 
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a 
b- 



1,2,3 
1 



(33) 



V 1 = {X',Wn - Uy(l) 

V 2 = (\ a ,WD -> SU L {2) 

% = {g\Wf) -> ST/ C (3) 

where PV/ 1 are the gauge fields and A', A and g are the and SUl{2) gauginos 

and the gluino, respectively. 

Leptons 

The leptons are described by chiral supermultiplets. As each chiral multiplet only 
describes one helicity state, we will need two chiral multiplets for each charged 
leptonP]. The multiplets are given in Table [|, where the t/y(l) hypercharge is defined 



Supermultiplet 


SU C (3) ® SU L (2) ® U Y (1) 
Quantum Numbers 


Li = (L, L)i 


(1,2,-1) 
(1,1,1) 



Table 3: Lepton Supermultiplets 

through Q = T s + Y. Notice that each helicity state corresponds to a complex scalar 
and that Li is a doublet of SUl(2), that is 



iiL 




(34) 



Quarks 

The quark supermultiplets are given in Table 
doublet of SUl(2), that is 



The supermultiplet Qi is also a 



Supermultiplet 


SU C (3) ® SU L {2) ® U Y (1) 
Quantum Numbers 


Qi = (Q, Q)i 

Di = (d R , d c L )i 
Ui = (u R ,u c L )i 


(3,2,|) 
(3,1,|) 
(3,1,-1) 



Table 4: Quark Supermultiplets 



Qi 




Qi 




(35) 



We will assume that the neutrinos do not have mass. 
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• Higgs Bosons 

Finally the Higgs sector. In the MSSM we need at least two Higgs doublets. This is 
in contrast with the SM where only one Higgs doublet is enough to give masses to 
all the particles. The reason can be explained in two ways. Either the need to cancel 
the anomalies, or the fact that, due to the analyticity of the superpotential, we have 
to have two Higgs doublets of opposite hypercharges to give masses to the up and 
down type of quarks. The two supermultiplets, with their quantum numbers, are 
given in Table [5| 



Supermultiplet 


SU C (3) ® SU L {2) ®U Y (1) 
Quantum Numbers 


Hi = (Hi, Hi) 
H 2 = (H 2 , H 2 ) 


(1,2,-|) 
(1,2,+D 



Table 5: Higgs Supermultiplets 



1.4.2 The Superpotential and SUSY Breaking Lagrangian 

The MSSM Lagrangian is specified by the R-parity conserving superpotential W 



W 



h^QtUjH b 2 + h%Q h fi 3 Hl + h%L\R 3 Ht - pHfH. 



(36) 



where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU(2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with £ 12 = 1. The coupling matrices hu^hn and He 
will give rise to the usual Yukawa interactions needed to give masses to the leptons and 
quarks. If it were not for the need to break SUSY, the number of parameters involved 
would be less than in the SM. This can be seen in Table ||. 

The most general SUSY soft breaking is 

V SB = MfQ a *Q a 3 - + MfUiU* + .u;f /;,/;; + Mfl^L] + MfttiR* + m 2 Hi H«*H« 
+m 2 H2 H^*H^ + [|M S A S A S + \M\\ + \M'X'X' + h.c. 



+£ab 



(37) 



1.4.3 Symmetry Breaking 

The electroweak symmetry is broken when the two Higgs doublets H\ and H 2 acquire 
VEVs 

El= (W + «+«S) t .„,) (38) 
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Theory 


Gauge 


Fermion 


Higgs 




Sector 


Sector 


Sector 


SM 


e,g,a s 


hu,ho, h-E 




MSSM 


e,g,a s 


hu,JlD, h-E 


n 


Broken MSSM 


e,g,a s 


hu,ho, h-E 


/i,M 1 ,M 2 ,M 3 ,A u ,A D ,A E ,B 

m H 2 , m H 1 , m Q, m U, m D, m h m R 



Table 6: Comparative counting of parameters 

with my? = \g 2 v 2 and v 2 = v\ + v\ = (246 GeV) 2 . The full scalar potential at tree level 
is 



dW 



dzi 



+ V D + V ao ft 



Vtotal — X] 

i 

The scalar potential contains linear terms 

V — +°rr° -4- T°rr° 

v linear — L l u l ' h 2 u 2 



(39) 



(40) 



where 



*i = ( m k +^)v 1 -Bfiv 2 + l{g 2 + g' 2 )v 1 {v 2 1 -v 2 2 ), 
h = im 2 H2 + n 2 )v 2 - Bfivi - \(g 2 + g' 2 )v 2 (v 2 - vj) 



(41) 



The minimum of the potential occurs for ti — (i — 1, 2). One can easily see that this 
occurs for m 2 H2 < 0. 

1.4.4 The Fermion Sector 

The charged gauginos mix with the charged higgsinos giving the so-called charginos. In 
a basis where ip +T = H 2 ) and ip~ T = the chargino mass terms in the 

Lagrangian are 

where the chargino mass matrix is given by 

M c = 



M 2 -4 



72^ 



T2 gV2 
H 



(43) 



and M is the SU(2) gaugino soft mass. The chargino mass matrix is diagonalized by two 
rotation matrices U and V defined by 



(44) 
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Then 

WMcV- 1 = Man (45) 
where Men is the diagonal chargino mass matrix. To determine U and V we note that 

Mq D = VM^McV 1 = U*M c Ml(U*) 1 (46) 

implying that V diagonalizes MqMc and U* diagonalizes McMq. In the previous expres- 
sions the Ff 1 are two component spinors. We construct the four component Dirac spinors 
out of the two component spinors with the conventions^, 

Xi = (Sri (47) 



In the basis i/j ot 
are given by 



—iX 3 , Hi, H$) the neutral fermions mass terms in the Lagrangian 



C m = --(r) T M N r + h.c. 



(48) 



where the neutralino mass matrix is 



M 



N 



Mi 







b' v 2 





M 2 


\gvi 


-\gv 2 




\gvx 







\g'v 2 


-\gv2 








(49) 



and Mi is the U(l) gaugino soft mass. This neutralino mass matrix is diagonalized by a 
4x4 rotation matrix N such that 



N* MjyN 1 = diag (m F o , m F o , m F o , m F o ) 



and 



(50) 
(51) 

The four component Majorana neutral fermions are obtained from the two component 
via the relation 

/ r"\ 

(52) 




1.4.5 The Higgs Sector 

In the MSSM there are charged and neutral Higgs bosons. Here we just discuss the neutral 



Higgs bosons. Some discussion on charged Higgs bosons is included in Section |3.2.2| . For 
a complete discussion see ref . || . In the neutral Higgs sector we have two complex scalars 
that correspond to four real neutral fields. If the parameters are real (CP is conserved in 
this sector) the real and imaginary parts do not mix and we get two CP-even and two 
CP-odd neutral scalars. The form of the mass matrices can be very much affected by the 
large radiative corrections due to top-stop loops and we will discuss both cases separately. 

4 Here we depart from the conventions of ref. [^) because we want the \~ to be the particle and not 
the anti-particle. 
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Tree Level 



The tree level mass matrices are 



M B 2 



cot/3 -1 ' 
— 1 tan (3 , 



±m|sin2/3 + 



' tan j3 — 1 
- 1 



and 



Mi 

2\ 



A 



where 



cot (5 J 
A = Bfi 



A 



(53) 
(54) 



/ tan (3 — 1 
V - 1 cot/3, 

Notice that <h t(Mj z ) = 0. In fact the eigenvalues of Mj 2 are and m\ = 2 A/ sin 2/3. The 
zero mass eigenstate is the Goldstone boson to be eaten by the Z°. A is the remaining 
pseudo-scalar. For the real part we have two physical states, h and H, with masses 



l 



"h,H — 2 

with the tree level relation 



m A + m z =F v { m \ + ^|) 2 — 4m 2 4 m| cos 2 2/3 



which implies 



2,2 2,2 

m/j + m H = m A + m z 



rrih < tua < mH 
m h < m z < ma 



(55) 
(56) 

(57) 



Radiative Corrections 



The mass relations in Eq. fl5?| ) were true before it was clear that the top mass is very 
large. The radiative corrections due to the top mass are in fact quite large and can 
not be neglected if we want to have a correct prediction. The whole picture is quite 
complicated ||, but here we just give the biggest correction due to top-stop loops. The 
mass matrices are now, in this approximation, 




sin 2/3 + 




(5f 



and 



where 



with 



A = Bfj,- 



M T 



3^ 2 




m: 



64tt 



2 sin (3 m 



f(m 2 



2m 2 



An 



i i 

VTLr — m~ 

*1 *2 




- 1 



(59) 



(60) 



(61) 
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The Aij are complicated expressions H]. The most important is 



mf , / ml m~ 



A 22 = log I (62) 

sin p \ I 



Due to the strong dependence on the top mass in Eq. Q6"2]) the CP-even states are the 
most affected. The mass of the lightest Higgs boson, h can now be as large as 140 GeU||. 



1.5 The Constrained Minimal Supersymmetric Standard Model 

We have seen in the previous section that the parameters of the MSSM can be consid- 
ered arbitrary at the weak scale. This is completely consistent. However the number of 
independent parameters in Table |G| can be reduced if we impose some further constraints. 
That is usually done by embedding the MSSM in a grand unified scenario. Different 
schemes are possible but in all of them some kind of unification is imposed at the GUT 
scale. Then we run the Renormalization Group (RG) equations down to the weak scale 
to get the values of the parameters at that scale. This is sometimes called the constrained 
MSSM model. 

Among the possible scenarios, the most popular is the MSSM coupled to N = 1 Super- 
gravity (SUGRA). Here at M GUT one usually takes the conditions: 

A t = A b = A T = A, B = A-l , 

m 2 Hl = m\ 2 = M 2 L = M\ = m\ ,M% = M 2 = M 2 D = m\ , 

M 3 = M 2 = Mi = Mi/2 (63) 

The counting of free parameters^ is done in Table [7]. 



Parameters 


Conditions 


Free Parameters 


ht, h b , hr, Vi, v 2 


m w , m t , m b , m T 


tan (3 


A, m , Mi/2, H 


t« = 0,z = l,2 


2 Extra free parameters 


Total = 9 


Total = 6 


Total = 3 



Table 7: Counting of free parameters in the MSSM coupled to N=l SUGRA 

It is remarkable that with so few parameters we can get the correct values for the param- 
eters, in particular m\ 2 < 0. For this to happen the top Yukawa coupling has to be large 
which we know to be true. 



J For one family and without counting the gauge couplings. 
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2 Spontaneous Breaking of R— parity 



2.1 Introduction 

Most studies of super symmetric phenomenology have been made in the framework of the 
MSSM which assumes the conservation of a discrete symmetry called R-parity (Rp) as 
has been explained in the previous Section. Under this symmetry all the standard model 
particles are R-even, while their superpartners are R-odd. R p is related to the spin (S), 
total lepton (L), and baryon (B) number according to R p = (—l)( 3B+L+2S \ Therefore the 
requirement of baryon and lepton number conservation implies the conservation of R p . 
Under this assumption the SUSY particles must be pair-produced, every SUSY particle 
decays into another SUSY particle and the lightest of them is absolutely stable. These 
three features underlie all the experimental searches for new supersymmetric states. 

However, neither gauge invariance nor SUSY require R p conservation. The most general 
supersymmetric extension of the standard model contains explicit R p violating interactions 
that are consistent with both gauge invariance and supersymmetry. Detailed analysis of 
the constraints on these models and their possible signals have been madeJF]. In general, 
there are too many independent couplings and some of these couplings have to be set to 
zero to avoid the proton to decay too fast. 

For these reasons we restrict, in this Section, our attention to the possibility that R p 
can be an exact symmetry of the Lagrangian, broken spontaneously through the Higgs 



mechanism |8|, |9|, |10[. This may occur via nonzero vacuum expectation values for scalar 
neutrinos, such as 

Vr = (vr t ) ; v L = (V Lt ) . (64) 

If spontaneous Rp violation occurs in absence of any additional gauge symmetry, it leads 
to the existence of a physical massless Nambu-Goldstone boson, called Majoron (J)f|. In 
these models there is a new decay mode for the Z° boson, Z° — > p + J, where p is a light 
scalar. This decay mode would increase the invisible Z° width by an amount equivalent to 
1/2 of a light neutrino family. The LEP measurement on the number of such neutrinos JTI] 



is enough to exclude any model where the Majoron is not mainly an isosinglet [|11|] . The 
simplest way to avoid this limit is to extend the MSSM, so that the R p breaking is driven 
by isosinglet VEVs, so that the Majoron is mainly a singlet 0. In this section we will 
describe in detail this model for Spontaneously Broken R-Parity (SBRP) and compare 
its predictions with the experimental results. 



2.2 A Viable Model for Spontaneous R— parity Breaking 

In order to set up our notation we recall the basic ingredients of the model for spontaneous 
violation of R parity and lepton number proposed in|5). The superpotential is given by 

W = h u QH u U + h d H d QD + h e LH d R 

+ (h H u H d -e 2 )<l> 

+ h u LH u v c + h$Sv c (65) 
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This superpotential conserves total lepton number and R p . The superfields ($, z/ c j, Si) are 
singlets under SU2 <8> U(l) and carry a conserved lepton number assigned as (0, —1, 1) re- 
spectively. All couplings h u , h d , h e , h u , h are described by arbitrary matrices in generation 
space which explicitly break flavor conservation. 

As we will show in the next section these singlets may drive the spontaneous violation of 
i? p ||, [12] leading to the existence of a Majoron given by the imaginary part of 

'v u H u - v d H d ) + —u T - —v c T + —S T (66) 



Vv 2X ' V V V 

where the isosinglet VEVs 

v R = (K) , v s = (S T ) (67) 



with V — Jv 2 R + Vg, characterize R p or lepton number breaking and the isodoublet VEVs 

v u = (H u ) , v d = (H d ) , v L = iy Lr ) (68) 
drive electroweak breaking and the fermion masses. 



2.3 Symmetry Breaking 
2.3.1 Tree Level Breaking 

First we are going to show that the scalar potential has vacuum solutions that break R p . 
Contrary to the case of the MSSM described in the previous section, the model described 
by Eq. (|65|) can achieve the breaking of SU(2) x U(l) at tree level, without the need of 
having some negative mass squared driven by some RG equation. The complete model 
has three generations and, as we will see, some mixing among generations is needed for 
consistency. But for the analysis of the scalar potential we are going to consider, for 
simplicity, a 1-generation model. 

Before we write down the scalar potential we need to specify the soft breaking terms. We 
write them in the form given in the spontaneously broken N = 1 supergravity models, 
that is 

V soft = m -Ah ®H u H d - Be 2 <& + Ch v v c vH u + Dh§v c S + h.c. 

+ m 2 u \H u \ 2 + m 2 d \H d \ 2 + m\\V\ 2 + m\\V c \ 2 + m 2 \S\ 2 + m 2 F \<£ \ 2 (69) 
At unification scale we have 

C=D= A ; B=A-2 

m 2 u = m 2 d = ■■■ = rh 2 (70) 

At low energy these relations will be modified by the renormalization group evolution. 
For simplicity we take C = D = A and B = A — 2 but let[] rp? u 7^ fh 2 d 7^ • ■ • 7^ rh^. Then 

6 Notice that for (H u ) ^ (Hd) we must have ^ fh d even in MSSM. 
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the neutral scalar potential is given by 
1 



V t 



total 



(g 2 + g 



I2\ 



H,.\ — \Hd 



\v\ 2 



+ \h&S + h v vH u ¥ + \h Q $H u r + \h<bi>' 



r.c 1 2 



+ \h$U' 



r.c 1 2 



- h ®H d + h v vv c \ l + | - h H u H d + hu c S - e 



2 1 2 



+ m [-A (-h$v c S + h $H u H d - h v vE u v c ) + (2 - A)s 2 <& + h.c. 



(71) 



where Zi stand for any of the neutral scalar fields. The stationary equations are then 

= . (72) 



dV 

dzi 



These are a set of six nonlinear equations that should be solved for the VEVs for each set 
of parameters. To understand the problems in solving these equations we just right down 
one of them, for instance 



dV 



dH rl 



H d =v d 



4 



+ 9 )K - v d - v l) - K v u -rh d - h v F 



(Ah m v F + hh v R v s - h e 2 )v u - h v v L v R h v F = (73) 



Also it is important to realize that it is not enough to find a solution of these equations 
but it is necessary to show that it is a minimum of the potential. To find the solutions 
we did not directly solve Eq. (|72|) but rather use the following three step procedure: 

1. Finding solutions of the extremum equations 

We start by taking random values for h, ho, h u , A, e 2 , m , vr, and i>g. Then choose 
tan (3 = v u /vd and fix v u , Vd by 



rn 



w 



12/2, 2 i 2 \ 

\9 K + v d + v L ) 



Finally we solve the extremum equations exactly for mj, m 2 



(74) 

mn. This is 



possible because they are linear equations on the mass squared terms. 

2. Showing that the solution is a minimum 

To show that the solution is a true minimum we calculate the squared mass matrices. 



These are 



Ml, 

1 IJ 



d 2 v , r \ , #y_ 



2 \ dzidzj 



i ( d 2 V | rr \ | dV 



(75) 
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The solution is a minimum if all nonzero eigenvalues are positive. A consistency 
check is that we should get two zero eigenvalues for M 2 corresponding to the Gold- 
stone boson of the Z° and to the majoron J. 

3. Comparing with other minima 

There are three kinds of minima to which we should compare our solution. 

• v u = v d = v L = v R = v s = ; v F ^0 

• v L = v R = v s = ; v u , v d , v F ^ 

• v u = v d = v L = ; v R , v s , v f ^ (76) 



As a final result we found a large region in parameter space where our solution that 
breaks R p and SU 2 ®U(1) is an absolute minimum. 



2.3.2 Radiative Breaking 

We tried to constrain the model of Eq. fl65|) by imposing boundary conditions at some 
unification scale and using the RG equations to evolve the parameters to the weak scale. 
Despite all our efforts we were not able to obtain radiative spontaneous breaking of both 
Gauge Symmetry and R p in this simplest model. 

To show the point that this could be achieved, we consider instead a model with Rank-4 
unification, given by the following superpotential: 

W = h u u c QH u + h d d c QH d + h e e c LH d 

+ h H u H d § + h u v c LH u + h$v c S + A$ 3 (77) 

The boundary conditions at unification are 

A u = A = A = A v = A x , 

Mk = M k = K = Ml = M 2 Q = ml, 

Ml = C v oml ; Ml = C s m 2 ; M| = C^m 2 , 

M 3 = M 2 = Mi = Mi/a (78) 

We run the RGE from the unification scale M\j ~ 10 16 GeV down to the weak scale. In 
doing this we randomly give values at the unification scale. After running the RGE we 
have a complete set of parameters, Yukawa couplings and soft-breaking masses m 2 (RGE) 
to study the minimization of the potential, 

Vtotai — 

i 

To solve the extremum equations we use the method described before: 



dW 

dzj 



+ V D + V SB + V RC 



(79) 
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1. The value of v u is determined from m top = h t v u for m top = 175 ± 5 GeV. If v u 
determined in this way is too high we go back to the RGE and choose another 
starting point. 

2. Vd and tan(/3) are then determined by mw 

3. vl is obtained by solving approximately the corresponding extremum equation. 

4. We then vary randomly m , Vr, v$, v^. 

5. We solve the extremum equations for the soft breaking masses, which we now call 

ml 

6. Calculate numerically the eigenvalues to make sure it is a minimum. 

After doing this we end up with a set of points for which: i) The Yukawa couplings and 
the gaugino mass terms are given by the RGE's, ii) For a given set of mf each point is 
also a solution of the minimization of the potential that breaks R p , in) However, the m\ 
obtained by minimizing the potential differ from those obtained from the RGE, m\[RGE) . 
Our goal is to find solutions that obey 

ml = m 2 t (RGE) Vi (80) 

To do that we define a function 



in: mj(RGE)\ 



V = max [ 2 'p . , '"' v '7 - ) V ^ ( 81 ) 
\mf(RGE) mf J 



From Eq. fl8"T| ) we can easily see that 77 > 1. We are then all set for a minimization 
procedure. We were not able to find solutions with strict universality. But if we relaxed[] 
the universality conditions on the squared masses of the singlet fields we got plenty of 
solutions. 



2.4 Main Features of the Model 

In this section we will review the main features of the model of spontaneous broken Rp 
described by Eqs. (EH) and (|69"D. 



2.4.1 Chargino Mass Matrix 

The form of the chargino mass matrix is common to a wide class of SUSY models with 
spontaneously broken R p and is given by |T^, [L3| 



7 This meant that the C's in Eq. ( |78|) were not equal to 1. A few percent of non-universality was 
enough to get solutions. 
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Hi 


-iW+ 






hvij ^Rj 




H d 




f 


gvd 


-iW~ 





gv u 


M 2 



(82) 



Two matrices U and V are needed to diagonalize the 5x5 (non-symmetric) chargino mass 
matrix 



xt = ; TP+ = (ef,e+,et,H+,-iW+) 

XT = UyipJ ; ipr = fc,e2,e£,Hj,-iW-) 
where the indices i and j run from 1 to 5. 



53) 



2.4.2 Neutralino Mass Matrix 

Under reasonable approximations, we can truncate the neutralino mass matrix so as to 



obtain an effective 7x7 matrix f 3 







H u 


H d 


-tW 3 


-iB 














g' 


H u 







-n 


q 

— t=V u 


I 


H d 





-fi 





q 

72 Vd 


g' 

-72 Vd 


-iW 3 




a 


q 


M 2 





-iB 


g' 

~V2 VLi 


g' 
V2 


g' 
V2 





Mi 
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This matrix is diagonalized by a 7 x 7 unitary matrix N, 

X° = where V° = K H u , H d , -iW 3 , -iB) 



2.4.3 Charged Current Couplings 



55) 



Using the diagonalization matrices we can write the charged current Lagrangian describing 
the weak interaction between charged lepton/chargino and neutrino/neutralinos as 



Cpc = J_ w ^- r{KL . k p L + K Rlk P R ) X l + h.c. 
where the 5x7 coupling matrices R may be written as 
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(86) 



Kuk = r]i(-\f2U i5 N k Q - U iA N kb - ^ U im N km ) 

m=l 

K Rlk = e k (-V2V l5 N k6 + V l4 N k4 ) (87) 



2.4.4 Neutral Current Couplings 



The corresponding neutral current Lagrangian may be written as 
r NC _ 9 

where the 7x7 coupling matrices 0' L R and 0" L R are given by 



COs6>i 



X7r(0' Lik P L + 0' Rlk P R ) Xk + \ X-nO'^PL + O' Rik P R )x k 



0' Lik = riirjl^Ui A U ki + U i5 U k5 + ^ H U im U km -5 ik sin 2 6w^J 
O'mk = 2 V i* V k4 + V i5 Vk 5 - S ik sin 2 9 W 

0"Lik = \ e i e k [Ni±N kA - iVjS^feS - X! N irnN km )=-€i€ k O Rik 
1 \ m=l / 



(89) 



In writing these couplings we have assumed CP conservation. Under this assumption the 
diagonalization matrices can be chosen to be real. The r]i and e k factors are sign factors, 
related with the relative CP parities of these fermions, that follow from the diagonalization 
of their mass matrices. 

2.4.5 Parameters values 

All the results discussed in the following sections use Eqs. (|86| ) and (|88| ) for the charged 
and neutral currents, respectively. To compare with the experiment we need to discuss 
the input parameters. Typical values for the SUSY parameters fj, = ho M 2 , and the 
parameters h U i^ lie in the range 

-250 < J^rr < 250 30 < t4% < 1000 

- GeV - - GeV - (9Q) 

io- 10 < h ul3 , h u23 < io- 1 io- 5 < h u33 < io- 1 

and we take the GUT relation M1/M2 = 5/3 tan 2 8w- For the expectation values we take 
the following range: 

v L = v L3 = 100 MeV v L1 = v L2 = 

50 GeV <v R = v R3 < 1000 GeV v R1 = v R2 = 

50 GeV <v s = v S3 = v R < 1000 GeV v si = v S2 = 

1 < tan B = — < 50 

v d 



(91) 
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which means that in practice we are considering that R p breaking is obtained only through 
r lepton number violation. 

2.4.6 Experimental Constraints 

Before we close this section on the spontaneously broken R p model we have to discuss 
what are the experimental constraints on the model. Some of these constraints are com- 
mon to all SUSY models, and are related to the negative results of the searches for the 
superpartners. This in turn puts constraints in the parameters of the models. But there 
are other constraints that are more characteristic of the spontaneously broken R p models, 
in particular those that are related to lepton flavor violation. We will give here a short 
list of the constraints that we have been using. 

• LEP searches 

The most recent limits on chargino masses from the recent runs were included. 

• Hadron Colliders 

From pp colliders there are restrictions on gluino production and hence on the gluino 
mass. 

• Non— Accelerator Experiments 

They follow from laboratory experiments related to neutrino physics, cosmology and 
astrophysics. The most relevant are: 

— Neutrinoless double beta decay 

— Neutrino oscillation searches 

— Direct searches for anomalous peaks at it and K meson decays 

— The limit on the tau neutrino mass 

— Cosmological limits on the v T lifetime and mass 

2.5 Implications for Neutrino Physics 

Here we briefly summarize the main results for neutrino physics. 

• Neutrinos have mass 

Neutrinos are massless at Lagrangian level but get mass from the mixing with 
neutralinos[llO|, |13 |. 



Neutrinos mix 

The coupling matrix h u . has to be non diagonal to allow 



Vt v„ + J (92) 
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and therefore evading |T3] the Critical Density Argument against v's in the MeV 
range. The fact that h Uij has to be non diagonal leads to important consequences 
in lepton violating processes as we will see below. 

Avoiding BBN constraints on the vh Vt 



In the SM BBN arguments [14] rule out v T masses in the range 



0.5 MeV < m Ur < 35 MeV (93) 



We have shown [13] that SBRP models can evade that constraint due to new anni- 
hilation channels 

z/ r z/ T -> J. J (94) 
2.6 R— parity in Non— Accelerator Experiments 

Here we will describe the implications of SBRP in non accelerator experiments like the 
solar neutrinos experiments and flavor violating leptonic decays. 

2.6.1 Solar Neutrinos 



To a good approximation we can write [13|] 



v\ = cos Qv e — sin Qv^ 
z/ 2 = sin Qv e + sin Qv^ 

z/ 3 = v r (95) 

where z/j, % = 1,2,3 and z/ e , and v T are, respectively, the mass and weak interaction 
eigenstates. The mixing angle 9 is given in terms of the model parameters by 

tan# = ^ (96) 

The constraints on h ui3 and h U23 do not restrict much their ratio. Therefore a large range 
of mixing angles is allowed. For the masses we get []13f 



mi = 
10~ 4 eV < m 2 < 10~ 2 eV 

lOkeV < m T < 23MeV (97) 



just in the right range for the MS W mechanism. 



2.6.2 SUSY Signals in jjt, and r Decays 



The existence of a massless scalar particle, the majoron, can affect the decay spectra of 
the \i and r leptons through the emission of the Majoron in processes such as 
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/i^e+J ; r^e + J ; r — > jj, + J (98) 

These are flavor violating decays that are present in our model because the matrix h V ij is 
not flavor diagonal. After a careful sampling of the parameter space, we found out that 
the rates can be close to the present experimental limits [|HJ . For instance for the process 
H — > eJ we can go up to the present experimental limit |T7[, BR(fi — > eJ) < 2.6 x 10 -6 . 

2.7 R-parity Violation at LEP I 
2.7.1 Higgs Physics 

The structure of the neutral Higgs sector is more complicated then in the MSSM. However 
the main points are simple. 

• Reduced Production 

Like in the MSSM the coupling of the Higgs to the Z° is reduced by a factor e# 



9ZZh 



9zzh 



< 1 (99) 



Invisible decay 

Unlike the SM and the MSSM where the Higgs decays mostly in bb, here it can have 
invisible decay modes like 

H^J + J (100) 

Depending on the parameters, the BR(H — > invisible) can be large. This will relax 
the mass limits obtained from LEP. We performed a model independent analysis of 



the LEP data (jTSj taking m#, Sb and BR(H invisible) as independent parame- 



ters. The results are shown in Fig. (|T]a) 

2.7.2 Chargino Production at the Z Peak 

The more important is the possibility of the decay 

z° - X ^ (101) 



This decay is possible because R p is broken. We have shown [|T0"| , |l9f that this branching 
ratio can be as high as 5 x 10~ 5 . This is shown in Fig. ([l|b). Another important point is 
that the chargino has different decay modes with respect to the MSSM. 

x - x° + fT 

X - r + J (102) 

The relative importance of the 2-body over the 3-body is very much dependent on the 
parameters of the model, but the 2-body can dominate. 
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2.7.3 Neutralino Production at the Z Peak 

We have developed an event generator that simulates the processes expected for the LEP 
collider at y/s = M z . Its main features are: 

• Production 

As far as the production is concerned, our generator simulates the following processes 
at the Z peak: 

e+e" -> x" (103) 



e + e" -> XX (104) 

Decay 

The second step of the generation is the decay of the lightest neutralino. The 2-body 
only contributes to the missing energy. The 3-body are: 

X -> v T Z* -> i/ T l + T,u T uu, v T q^ (105) 



X -> r TV* -> TUik, rq u q d (106) 

• Hadronization 

The last step of our simulation is made calling the PYTHIA software for the final 
states with quarks. 
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One of the cleanest and most interesting signals that can be studied is the process with 
missing transverse momentum + acoplanar muons pairs [2(J 



h + /iV (107) 

The main source of background for this signal is the 

Z — > ji + n~ + soft photons (108) 

For defmiteness we have imposed the cuts used by the OPAL experiment for their search 
for acoplanar dilepton events: (a) We select events with two muons with at least for 
one of the muons obeying | cos^l less than 0.7. (b) The energy of each muon has to be 
greater than a 6% of the beam energy, (c) The missing transverse momentum in the event 
must exceed 6% of the beam energy, j> T > 3 GeV. (d) The acoplanarity angle (the angle 
between the projected momenta of the two muons in the plane orthogonal to the beam 
direction) must exceed 20°. With these cuts we were able to calculate the efficiencies of 
our processes. 

We used the data published by ALEPH in 95 and analyzed both the single production 
e + e~ — ► xv and the double production e + e~ —> XX processes. For single production we 
get 

N exp t(xv) = <y{e + e~ -> x^)BR(x -> v T ^ + ^i~)e xu L int (109) 

Using the expression for the cross section we can write this expression in terms of the 
product BR{Z — ► x v ) x BR(x v T fi + fi~) and obtain a 95%CL limit on this R-parity 
breaking observable, as a function of the x mass. This is shown in Figure £| For the 
double production of neutralinos the number of expected j>T + events is 

N exp t(xx) = cr(e + e~ -> xx)2BR(x -> invisible) BR(x -> v T iJ + n~)e xx L int (110) 

We can obtain an illustrative 95%CL limit on BR(Z — > xx) x BR{x ~ * v T [i + n~) x 
BR(x —>■ invisible) as a function of the x mass P0| . This is also shown in Figure |2] where 
we can see that the models begin to be constrained by the LEP results. 



2.8 R-parity Violation at LEP II 
2.8.1 Invisible Higgs 



The previous LEP I analysis has been extended for LEP II.|2T| As a general framework 
we consider models with the interactions 

Chzz = e B (V2G F ) 1/2 M 2 z Z^h, 

Chaz = -tA—^—Z^hd^A, (111) 
cos 6w 

with €a(b) being determined once a model is chosen. We also consider the possibility that 
the Higgs decays invisible 

h^JJ (112) 
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Figure 2: On the left a comparison of the attainable limits on BR(Z — * x v )BR(x — > 
versus the lightest neutralino mass, with the maximum theoretical values expected in different R 
parity breaking models. The solid line (a) is just for the /i + '\x~v channel, while (b) corresponds to 
the improvement expected from including the e + e~v channel, as well as the combined statistics 
of the four LEP experiments. The dashed line corresponds to a model with explicit R-parity 
violation, while the dotted one is calculated in the spontaneous R-parity-violation model. On 
the right the same for BR(Z — > xx)BR{x — * u + [i~v). 



e e 



e e 



(113) 



and treat the branching fraction B for h — > J J as a free parameter. 
The following signals with p T were considered: 

-> (Zh + Ah) ^bb + f T , 

_, zh^ £+r + f T , 

but also the more standard processes 

e + e~ -> Zh^£ + £~+ 66, 

e+e- -> (Z/i + A/i) ^bb + 66 . (114) 
Using the above processes and after a careful study of the backgrounds and of the necessary 



cuts, [f2T[ it was possible to evaluate the limits on M^, e^, es, and B that can be 
obtained at LEP II. In Figure |3] are shown some of these limits. 



2.8.2 Neutralinos and Charginos 

At LEP II the production rates for R-parity violation processes will not be very large, 
compared with those at LEP I. Therefore we expect that the production rates will be like 
in the MSSM, via non R-parity breaking processes. However the decays will be modified 
much in the same way as in the LEP I case. This is specially important for the xo because 
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Figure 3: 
bounds on 



On the left, bounds on e 2 B as a function of Mh for \fs = 175GeV. On the right, 
e\ as a function of Mh and Ma for B = 1 and y/s = 175 GeV. 



it is invisible in the MSSM but visible here. Also the R-parity violating decays of the 
charginos 

X'^r- + J (115) 
can have a substantial decay fraction compared with the usual MSSM decays 

x-^x° + fT (116) 
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3 Bilinear R— parity Violation: The e model 



We have seen in the previous section that it could well be that R-parity is a symmetry 
at the Lagrangian level but is broken by the ground state. Such scenarios provide a 
very systematic way to include R parity violating effects, automatically consistent with 
low energy baryon number conservation. They have many added virtues, such as the 
possibility of providing a dynamical origin for the breaking of R-parity, through radiative 
corrections, similar to the electroweak symmetry p2| . The simplest truncated version of 
such a model, in which the violation of R-parity is effectively parameterized by a bilinear 
superpotential term eiL\H\ has been widely discussed |23], It has also been shown 
recently |24[] that this model is consistent with minimal N=l supergravity unification with 
radiative breaking of the electroweak symmetry and universal scalar and gaugino masses. 
This one-parameter extension of the MSSM-SUGRA model therefore provides the simplest 
reference model for the breaking of R-parity and constitutes a consistent truncation of 
the complete dynamical models with spontaneous R-parity breaking proposed previously 
0. In this case there is no physical Goldstone boson, the Majoron, associated to the 
spontaneous breaking of R-parity, since in this effective truncated model the superfield 
content is exactly the standard one of the MSSM. Formulated as an effective theory at 
the weak scale, the model contains only two new parameters in addition to those of 
the MSSM. Therefore our model provides also the simplest parameterization of R-parity 
breaking effects. In contrast to models with tri-linear R-parity breaking couplings, it leads 
to a very restrictive and systematic pattern of R-parity violating interactions, which can 
be taken as a reference model. In this section we will review the most important features 
of this model. 



3.1 Description of the Model 

The superpotential W is given by 



W 



hZQtUjH* + hZQ'DjH? + ti^RjHl - ^H b + e t L*H b (117) 



where i, j = 1, 2, 3 are generation indices, a, b = 1, 2 are SU(2) indices. In the following we 
will consider, for simplicity, only the third generation. Then the set of soft supersymmetry 
breaking terms are 

V sof t = M 2 Q Ql*Ql + M 2 U£U 3 + M 2 D D;D 3 + MllfLl + M 2 R R* 3 R 3 + m 2 Hl H?H; 



+ m 2 H2 H^m- 



±M 3 A 3 A 3 + |M 2 A 2 A 2 + 5M1A1A1 + h.c. 



+ tab 



A t h t Q a 3 U 3 H b +A b h b Q b 3 D 3 H? +A T h T L b 3 R 3 H*-Bu,H«H b +B 2 e 3 L a 3 H b 



-a zrb 



118) 



The bilinear R p violating term cannot be eliminated by superfield redefinition. The reason 
is that the bottom Yukawa coupling, usually neglected, plays a crucial role in splitting the 
soft-breaking parameters B and B2 as well as the scalar masses m 2 H and M£, assumed to 
be equal at the unification scale. 
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The electroweak symmetry is broken when the VEVS of the two Higgs doublets Hi and 
if 2, and the tau-sneutrino. 



Hi 



V 



V2 



, H 2 = 



( Hi \ 



( v? + v 3 + iK \ 



r.I 



J 



V 



V2 



) 



V 



V2 



(119) 
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The gauge bosons W and Z acquire masses = \g 2 v 2 , m 2 z = |(g 2 + g' 2 )v 2 , where 

v 2 = v\ + v\ + v\ = (246 GeV) 2 (120) 
We introduce the following notation in spherical coordinates: 



v i = v sin 9 cos (5 
v 2 — v sin 9 sin (5 

V 3 = V COS 9 



(121) 



which preserves the MSSM definition tan/3 = v 2 /vi. The angle 9 equal to tt/2 in the 
MSSM limit. 



The full scalar potential may be written as 



Vtotal — X/ 



dW 



+ V D + Vsoft + V RC 



(122) 



where Zi denotes any one of the scalar fields in the theory, Vd are the usual D-terms, 
V so f t the SUSY soft breaking terms, and Vrc are the one-loop radiative corrections. In 
writing Vrc we use the diagrammatic method and find the minimization conditions by 
correcting to one-loop the tadpole equations. This method has advantages with respect 
to the effective potential when we calculate the one-loop corrected scalar masses. The 
scalar potential contains linear terms 



T/ _ M U , M U , J.0 -R 

^linear — L IA1 ' L 2A2 ' l 3 u t ' 

(m 2 Hi + /i 2 )^i - B/j,v 2 - /J,e 3 v 3 + \(g 2 + g' 2 )vi(vf -v 2 + v\) , 
(m 2 H2 + fi 2 + e\)v 2 - BfjLVt + B 2 e 3 v 3 - \{g 2 + g' 2 )v 2 (v 2 - v\ + v 2 3/ 
( m l 3 + 4) v z ~ A*e 3 ui + B 2 e 3 v 2 + \(g 2 + g' 2 )v 3 {v{ -v\ + v\) . 



where 



'i 



(123) 



(124) 



These t®,i = 1, 2, 3 are the tree level tadpoles, and are equal to zero at the minimum of 
the potential. 
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3.2 Main Features 



3.2.1 Charginos and Neutralinos 

The e-model is a one parameter generalization of the MSSM. It can be thought as an 
effective model showing the more important features of the SBRP-model at the weak 
scale. In fact the mass matrices, the charged and neutral currents, are similar to the 
SBRP-model if we identify 

£3 = V R hu33 (125) 

Therefore all that we said about the SBRP-model in Section 2 also applies here. In 
particular the implications of the mixing of the r lepton with charginos have been studied 
in ref. [25]. Their results are shown in Fig. |]a, and are similar to those of Section 2.7.2 if 



we use the identification of Eq. ( |125| ). The only difference arises in processes where the 
Majoron plays an important role, because it is absent here. This has been studied in full 
detail in refs. [E0|, 121 



The other important feature it is that this model has the MSSM as a limit. This can be 
illustrated in Fig. |4]b, where we show the ratio of the lightest CP-even Higgs boson mass 
rrih in the e-model and in the MSSM as a function of v^. As v 3 — > the ratio goes to one. 




65 70 75 80 85 90 95 20 40 60 80 100 



raft (GeV) |v 3 | (GeV) 

a) b) 

Figure 4: a) Regions of attainable cross section in BRpV in the plane tau neutrino mass vs 
chargino mass including large values of tan /3. b)Ratio of the lightest CP-even Higgs boson mass 
in the e-model and in the MSSM function of v 3 . 



3.2.2 Charged Scalars 

The charged scalar sector is also similar to the SBRP-model, because the extra superfields 
needed in that case are all neutral. Therefore the charged scalars are the charged Higgs 
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bosons, sleptons and squarks. Because of the breaking of R-parity the charged Higgs 
bosons are mixed with the charged sleptons. If we consider only the third generation, the 
mixing will be with the staus. Although this sector is similar in the e-model and in the 
SBRP-model, the overall analysis is simpler in e-model because it has fewer parameters. 

The mass matrix of the charged scalar sector follows from the quadratic terms in the 
scalar potential 



V, 



quadratic 



[Hi ,H2,f£,fn]M^ 



T 



T 



R 



+ 



(126) 



For convenience reasons we will divide this 4x4 matrix into 2x2 blocks in the following 

way: 



M 2 HH 



M 



2 

ff 



(127) 



where the charged Higgs block is 



M 2 HH - 



B l if + \g\vl-vl)+^ + \hlvl + ^ 



Bfi + \g 2 v x v 2 



Bp+\?W2 B^ + \g 2 (v 2 + v 2 )-B 2 e 3 f 2 +^ 



;i28) 

and h T is the tau Yukawa coupling. This matrix reduces to the usual charged Higgs mass 
matrix in the MSSM when we set v 3 = e 3 = and we call m\ 2 = B\x. The stau block is 
given by 



1l2„,2 1^2/ 



SKvt 



B 2 e 3 ^ + ^ 



-^h T (A T vi - pbv 2 ) 



(129) 

We recover the usual stau mass matrix again by replacing t> 3 = e 3 = 0, nevertheless, we 
need to replace the expression of the third tadpole in Eq. (|124 ) before taking the limit. 
The mixing between the charged Higgs sector and the stau sector is given by the following 
2x2 block: 



M 



Hf 



\h\v ii> 3 + lg 2 V!V 3 



- -j^h T (e 3 v 2 + A T 



v 3j 



-B 2 e 3 + \g 2 v 2 v 3 
-^h T (fiv 3 + e 3 vx) 



(130) 



and as expected, this mixing vanishes in the limit v 3 = e 3 = 0. The charged scalar mass 
matrix in Eq. (|127 ), after setting t\ — t 2 — t 3 — 0, has determinant equal to zero since 
one of the eigenvectors corresponds to the charged Goldstone boson with zero eigenvalue. 

The numerical study of the lowest-lying charged scalar boson mass has been done in 
ref. p6 |. The results are illustrated in Fig. [5|a. The main point to note is that can be 
lower than expected in the MSSM, even before including radiative corrections. This is due 
to negative contributions arising from the R-parity violating stau-Higgs mixing, controlled 
by the parameter e 3 . An alternative way to display the influence of e 3 parameter on the 
charged Higgs boson mass can be seen in Fig. ^b. In this figure the curves corresponding 
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25 50 75 100 125 150 0.5 1 5 10 50 100 



m A (GeV) tan/S 

a) b) 

Figure 5: a) Tree level and one-loop charged Higgs boson mass as a function of the CP-odd 
Higgs mass tua- The horizontal dashed line corresponds to the TV-boson mass, b) Minimum 
of the charged Higgs boson mass versus tan (5. Each curve corresponds to a different range of 
variation of the R-parity violating parameters £3 and 113. 



to different £3 and t>3 values delimit the minimum theoretically allowed charged Higgs 
boson mass corresponding to those specific values. 

We now turn to a discussion of the charged scalar boson decays. In Fig. |]a we display [^6 



the stau decay branching ratios below and past the neutralino threshold and in Fig. |6|b 
the charged Higgs branching ratios possible in the model for a particular set of chosen 
parameters. Finally, for the case of the R-parity violating charged Higgs boson decays 
one can see from Fig. |6|b that the branching ratios into supersymmetric channels can be 
comparable or even bigger than the R-parity conserving ones, even for relatively small 
values of e and V3. Another way to see that the dominance of R-parity- violating Higgs 
boson decays is not an accident of the above parameter choice is illustrated in Fig. [7]. The 
various curves denote the maximum attainable values for the R-parity-violating Higgs 
boson branching ratio B(H + — > t + Xi)- 



3.3 Radiative Breaking 

In the previous discussion of the e-model the parameters were varied at the weak scale 
with no restrictions besides the experimental constraints on the masses of the particles. 
However, as we have seen with the MSSM, the parameter space can be constrained if we 
embed the theory in a grand unified scenario. This can also be done in the e-model, both 



with [23] and without [2^] b — r unification. We will describe below these two possibilities. 



3.3.1 Radiative Breaking in the e model: The minimal case 

At Q = M GUT we assume the standard minimal supergravity unifications assumptions, 
A t = A b = A T = A ; B = B 2 = A - 1 , 
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Figure 6: a) Stau branching ratios possible in our model for a particular choice of parame- 
ters. Note the neutralino threshold below which only R-parity violating decays are present, b) 
Charged Higgs branching ratios possible in our model for a particular choice of parameters. 
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In order to determine the values of the Yukawa couplings and of the soft breaking 
masses at low energies we first run the RGE's from the unification scale Mqut 
GeV down to the weak scale. We randomly give values at the unification scale 
parameters of the theory. 
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(132) 



The value of h%. GUT /An is defined in such a way that we get the r mass correctly. As the 
charginos mix with the tau lepton, through a mass matrix is given by 

M -±gv u ' 



Mr 



T2 gVd 
1 



V2 



gv 3 



-e 3 



y/2' 
j= 2 h T v d 



Imposing that one of the eigenvalues reproduces the observed tau mass m T , h T 
solved exactly as [23| 



h 



2m 2 T 

Vd 



1 + Si 
l + 5 2 



(133) 

can be 
(134) 
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Figure 7: The curves denote the maximum attainable i? p -violating charged Higgs branching 
ratio versus tan j3. 



where the Si, i = 1,2, depend on m T , on the SUSY parameters M, /i, tan f3 and on the R p 



violating parameters £3 and V3. It can be shown that [24 



lim 5j 





«3 



(135) 



After running the RGE we have a complete set of parameters, Yukawa couplings and soft- 
breaking masses mf(RGE) to study the minimization. This is done by using a method 
similar to the one described before in Section 2: 



1. We start with random values for h t and at Mqut- The value of h T at Mqut is 
fixed in order to get the correct r mass. 

2. The value of V\ is determined from rrib = hbVi/y/2 for = 2.8 GeV (running b 
mass at mz)- 

3. The value of t>2 is determined from m t = h t V2/V2 for m t = 176 ± 5 GeV. If 

4 

vj + v\ > v 2 = — mL = (246 GeV) 2 (136) 

g 

we go back and choose another starting point. 

4. The value of V3 is then obtained from 



We see that the freedom in h t and hf, at Mqut can be translated into the freedom in 
the mixing angles f3 and 9. Comparing, at this point, with the MSSM we have one extra 
parameter 6. We will discuss this in more detail below. In the MSSM we would have 
9 = tt/2. 
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After doing this, for each point in parameter space, we solve the extremum equations, 
for the soft breaking masses, which we now call mj (i = Hi, H2, L). Then we calculate 
numerically the eigenvalues for the real and imaginary part of the neutral scalar mass- 
squared matrix. If they are all positive, except for the Goldstone boson, the point is a 
good one. If not, we go back to the next random value. As before, we end up with a set 
of solutions for which the mf obtained from the minimization of the potential differ from 
those obtained from the RGE, which we call mf(RGE). Our goal is to find solutions that 
obey 

mi = (RGE) Vz (138) 

To do that we define a function 



2 2 1 

mf mf 



"= ma H^W'^F) Vt (139) 

that satisfies 77 > 1. Then we are all set for a minimization program. For this we used 
the CERN Library Program MINUIT. Following this procedure we were able to find [^4j 
plenty of solutions. 

Let us discuss the counting of free parameters in this model and in the minimal N=l 
supergravity unified version of the MSSM. In the MSSM we have the parameters shown 
in Table [5]. Normally the two extra parameters are taken to be the masses of the Higgs 
bosons h and A, the lightest CP-even and the CP-odd states, respectively. For the e- 
model the situation is described in Table |8|. As we have said before there is an extra 
parameter. Finally, we note that in either case, the sign of the mixing parameter fi is 
physical and has to be taken into account. 



Parameters 


Conditions 


Free Parameters 


h t , h b , h T , v x , v 2 , v 3 


m w , m t , m b , m T 


tan/5, cos 9 


A, m , Mi/2, H, e 3 


ti = 0, % = 1,2,3 


2 Extra free param. 


Total =11 


Total = 7 


Total = 4 



Table 8: Counting of free parameters in the e-model 



3.3.2 Gauge and Yukawa Unification in the e model 

Besides achieving gauge coupling unification, GUT theories also reduce the number of 
free parameters in the Yukawa sector. In SU(5) models, h b = h T at M GUT . The predicted 
ratio m b /m T at My/EAK agrees with the experimental values. In the MSSM a relation 
between m top and tan (3 is predicted. Two solutions are possible: low and high tan (3 . In 
SO (10) and E% models h t = h b = h T at Mqut- In this case, only the large tan/3 solution 
survives. Recent global fits of low energy data (B(b — > 57) and the lightest Higgs mass) to 
the MSSM show that it is hard to reconcile these constraints with the large tan (3 solution. 
Also the low tan f3 solution with \i < is disfavored. 
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Motivated by these considerations we analyzed the gauge and Yukawa unification in the 
e-model. We found ]27| that the e-model allows b — r Yukawa unification for any value of 



tan /3 and satisfying perturbativity of the couplings. We also found the t — b — r Yukawa 
unification easier to achieve than in the MSSM, occurring in a wider high tan/3 region. 
We will describe below how we got these results. 

As before h T can be solved exactly 



2m 1 

Vd 



l + Si 
l + S 2 



(140) 



where the 5j , i — 1, 2, depend on m T , on the SUSY parameters M, /z, tan f3 and on the R p 
violating parameters e% and V3. Also h t and h b are related to m t and m b 

v v 
mt = h t — y= sin (3 sin 9 , = hf,— y= cos/3 sin 9 (141) 

where 

v = 2mw/g tan/3 = t> M /t>d cos9 = v 3 /v (142) 

In our approach we divide the evolution into three ranges: i) From mz — > m t we use 
running fermion masses and gauge couplings, ii) From m t — > Msusy w e use the two-loop 
SM RGE's including the quartic Higgs coupling A. in) Finally from M S usy Mqut 
we use the two-loop RGE's. Using a top — > bottom approach we randomly vary the 
unification scale Mqut and the unified coupling acur looking for solutions compatible 
with the low energy data 

a-^{m z ) = 128.896 ±0.090 
sin 2 9 w (m z ) = 0.2322 ± 0.0010 

a s (m z ) = 0.118 ±0.003 (143) 

We get a region centered around M GUT w 2.3 x 10 16 GeV acur 1 ~ 24.5 Next we use 
a bottom — > top approach to study the unification of Yukawa couplings using two-loop 
RGEs. We take 

m w = 80.41 ± 0.09 GeV 
m T = 1777.0 ± 0.3 MeV 

m b (m b ) = 4.1 to 4.5 GeV (144) 

We calculate the running masses m T (m t ) = r]~ 1 m T (m T ) and mj(mt) = ?7 b " 1 mfc(m{,) where 
t] T and r)b include three-loop order QCD and one-loop order QED. At the scale Q = m t 
we keep as a free parameter the running top quark mass m t (m t ) and vary randomly 
the SM quartic Higgs coupling A. In doing the running we used the following boundary 
conditions: 

1. At scale Q = m t 

A 2 (m t ) = 2m 2 i {m t )/v 2 ; i = t,b,r (145) 
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2. At scale Q = M SUS y 



M m susy) = h t {M^ USY ) sin (3 sin 9 
^b(Mg USY ) = h b (M£ USY ) cos sin 9 

^A m susy) = h A M susy) cos/? sin e Jj^J^ 



(146) 



where hi denote the Yukawa couplings of our model and Aj those of the SM. The boundary 
condition for the quartic Higgs coupling is 



A(M 



SUSY J 



(9 2 (M+ USY ) + g' 2 (M+ USY )] (cos 2/3 sin 2 9 + cos 



(147) 



The MSSM limit is obtained setting 9 — > tt/2 i.e. = 0. 

The results are summarized in Fig. H. The dependence of our results on a s and m b is 
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Figure 8: Top quark mass as a function of tan/3 for different values of the R-Parity violating 
parameter V3. Bottom quark and tau lepton Yukawa couplings are unified at Mqut- The 
horizontal lines correspond to the la experimental mt determination. Points with t — b — r 
unification lie in the diagonal band at high tan/3 values. We have taken M$usy = m t- 

totally analogous to what happens in the MSSM. The upper bound on tan/3, which is 
tan/3 < 61 for a s = 0.118, increases with a s and becomes tan/3 < 63 (59) for a s = 0.122 
(0.114). The top mass value for which unification is achieved for any tan/3 value within 
the perturbative region increases with a s , as in the MSSM. As for the dependence on 
m b , if we consider rrib(m b ) = 4.1 (4.5) GeV then the upper bound of this parameter is 
given by tan/3 < 64 (58). In addition, the MSSM region is narrower (wider) at high tan/3 
compared with the mb{mb) = 4.3 GeV case. The line at high tan/3 values corresponds to 
points where t — b — r unification is achieved. Since the region with [-^3 1 < 5 GeV overlaps 
with the MSSM region, it follows that t — b — r unification is possible in this model for 
values of \v 3 \ up to about 5 GeV, instead of 50 GeV or so, which holds in the case of 
bottom-tau unification. 
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